Covariant description of spin 1/2 particles in self— action 
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The bispinor wave function finds its fundamental application in the study of electrons, neutrinos 
and protons as particles bound by their own potentials. Classical electromagnetism and the Dirac 
electron theory appear to be natural extensions of the physics that rules the electron itself. However, 
the charge density and the negative energy eigenvalues are two concepts which cannot survive further 
scrutiny. The electron and the proton resemble minute planetary systems where the bispinor terms 
defining their physical properties revolve with definite angular momentum around the singular point 
of the particles' potentials. The mass of the electron-neutrino in relation to the electron's and the 
mass of the latter in relation to the proton's are determined under particular assumptions. The 
most important feature of the new approach is the faithful representation of free particles at rest 
with respect to inertial systems of reference. This certainly enhances the foundation of physics; 
precisely in the realm of Quantum Mechanics, where the customary description of free electrons is 
a confusing consequence of the limited scope of the Schrodinger and Dirac equations. 

PACS numbers: ll.10.Cd, 14.60.Cd 



I. GENERAL CONSIDERATIONS 

The free particle solutions of the Schrodinger equation, 
plane waves extending throughout space, bear no direct 
relation with the concept of particle that most of us may 
have. It is only natural that this circumstance gave rise 
to philosophical currents and crosscurrents from which 
no real consensus was ever obtained: on examining the 
probabilistic interpretation of the wave function, Niels 
Bohr arrived at the conclusion that science not necessar- 
ily should be concerned on whether or not the dynamic 
variables of the electron coexist with definite values. In 
essence, he argued that physical reality cannot be disso- 
ciated from the process of measurement. Consequently 
free electrons, which by definition do not interact even 
with measuring devises, remain metaphysical until their 
state of isolation is disrupted by experiment. Moreover, 
the Heisenberg uncertainty relations have revealed an in- 
trinsic limitation in the accuracy that can be attained 
when the position and the momentum of the electron are 
measured simultaneously. Hence the entire formalism of 
Quantum Mechanics is to be considered as the tool for 
deriving predictions of statistical nature Q 

However, contrary to the familiar answer of the Copen- 
hagen school, perhaps in accordance with Einstein's faith 
in an objective reality, physics is sustained on the princi- 
ple of inertia, which does not allow room for indetermi- 
nacies of the kind. For instance, the parameter of relative 
velocity appearing in the Lorentz transformations implies 
the presence of a particle at rest with respect to one of the 
inertial systems, and in uniform motion with respect to 
the other. Why? Because the relative velocity between 
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systems of coordinates is devoid of physical meaning un- 
til a real object (particle) appears in space. We may ask, 
what sense does it make to verify the Lorentz covariance 
of the Dirac equation when the equation itself cannot de- 
scribe a particle in isolation in accordance with the prin- 
ciple that makes the Lorentz transformation meaningful: 
the Galileo-Newton principle of inertia in the first place. 

As a matter of fact, Erwin Schrodinger solved the hy- 
drogen atom thinking of a proton fixed at the origin of an 
inertial system implied by the potential r _1 in the equa- 
tion. If the electron and the proton revolve around the 
center of mass, it is irrelevant for the point in question: 
it only means that in quantum mechanics one assumes 
that the atom remains fixed with respect to an inertial 
system. It follows, therefore, that in the attempt to pre- 
serve its conceptual consistency, quantum mechanics had 
to invoke the process of measurement or to rely on some 
other interpretation. Physics evolved leaving behind im- 
portant conceptual problems unsolved. And so, some 
twenty years after the advent of the Schrodinger equa- 
tion this problem emerged with a different face. Physi- 
cists became fully aware that, in order to complete the 
theory of the free electron, self-action had to be reck- 
oned with. Then QED offered an approach based on the 
idea that the point-electron emits and absorbs virtual 
and real photons along its path. This time, however, the 
divergent self-energy put a definite end to the hope of a 
consistent theory, not to mention that the QED approach 
itself was inconsistent with the principle of inertia from 
the very outset, since the free electron path cannot be 
made to vanish with the choice of system of reference • 

Free electrons exist independently of experiment; we 
think that it is of the concern of physics to explain their 
existence in a way consistent with general principles. For- 
tunately, the framework to analyze the problem is defined 
already: the description of a self-interacting electron at 
rest with respect to an inertial system of reference should 
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be carried out with the solutions of the one and only fun- 
damental equation that can be proposed for the study of 
spin 1/2 particles in self-action: 



+ Su,(e,u,p)} -0 = 0, 



(1) 



where Sn represents the self-action 4-vector of the par- 
ticle under study. 

The conceptual structures of equation (1) and of the 
Dirac equation 



In [dju + aAfj] tp = -i(mch x ) 7 5 ip, 



(2) 



may reinforce each other in their own range of applica- 
bility and may fuse at some point too. The free particle 
solutions of the customary wave equations are important, 
but only as part in the solution of problems with discon- 
tinuous external interaction, either in space or in time. 
The tunnel effect -through the potential barriers that 
limit the force free region of the potential well- is a typ- 
ical case where plane waves satisfying boundary condi- 
tions become necessary and meaningful also. Although, 
we suspect that plane waves with no boundary conditions 
are meaningless solutions which have misled physicists for 
over 75 years now. 



II. THE ELECTRON SELF-POTENTIALS 

Mostly every physicist knows that the energy of a point 
electron diverges; although not every physicist is well 
aware that Gauss' law cannot hold when the particle is 
deprived of physical dimensions. Indeed, the volume in- 
tegral of the last two expressions below reveals a contra- 
diction: 

V E = r~ 2 d r (r 2 r • E) = r~ 2 d r (l) = 4?r5(r). (3) 

The reason some textbooks invoke the divergence theo- 
rem to prove © correct is easy to understand: empty 
space cannot be the origin of the Coulomb potential. In- 
teresting enough, this is not the only case where physi- 
cal problems have been circumvented at the expense of 
mathematical congruence. Paul Dirac himself had the 
opinion that renormalization has defied all the attempts 
of the mathematician to make it sound Q; yet, such pro- 
cedures have advocates and detractors who may never 
convince one another. 

It is appropriate to recall that the functional depen- 
dence of the potential on the charge density brings about 
problems of its own. For example, who can tell if the 
electrostatic energy of an electron is localized where the 
potential and the charge density coexist or is it dispersed 
throughout space? 



£ = 



p4> 



V4>-V(f>. 



(4) 



More important, the association of the square of the elec- 
trostatic field with the energy density leads to the follow- 
ing conclusion: the ratio of electromagnetic momentum 



over electromagnetic energy of an electron in uniform mo- 
tion is incompatible with relativity p|. 

Experience shows that each and every electron in 
a given distribution contributes with their individual 
Coulomb potential to the total potential. However, if 
the electron happens to be a particle with physical di- 
mension, there is no reason to think that a small portion 
of its volume would contribute with a small portion to 
its Coulomb potential. Again, since there is no consistent 
model of the electron, the physics that make its existence 
possible is not known. 

The vanishing of the Laplacian of the inverse distance 
occurs in three dimensions only: 
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(5) 



In fact, this is a mathematical characteristic of the 
Coulomb potential inherent to the electron. This is a pos- 
tulate and it is incompatible with the concept of charge 
density. 

The task now is to adjust the electromagnetic theory 
to a new conceptual structure. We shall do that in parts 
and as the circumstance allows. It is convenient to say in 
advance that suitable combinations of the products of the 
four components of the wave function and their complex 
conjugates enable the construction of covariant densities. 
These densities play an essential role in the definitions 
of the electromagnetic inertia and of the other proper- 
ties of the electron; all of which properties will vanish 
beyond a small radius from the singular point of the po- 
tential (from now on "singularity" for short). In this way 
the Coulomb potential springs up from the inside of the 
electron; though, no functional relation with any of the 
covariant densities exists. 

Electromagnetic signals from the singularity necessar- 
ily relate to the 4-vector of null length, r a = (r, 7-) . The 
contraction of r a with the 4- velocity of the singular point, 

u a = (GUMl-u 2 )- 1 ^ -> (0,1), (6) 

yields the fundamental invariant involved in the Lienard- 
Wiechert potentials 

Io = -(raUa)" 1 -> r _1 , (7) 
A a = l u a = (u, 1) (r - ?• u)- 1 -> (O,^ 1 ), 

(8) 

(the arrow reads -for a particle at rest reduces to-). 

Direct consequence of their genesis, these potentials 
satisfy the fundamental differential equations 



d a d a Af} = 0, 
d a A a = 0, 



(9) 
(10) 



Consequently, the field tensor F a p — d a Ap — dp A a sat- 
isfies Maxwell's equations with no sources: 

OaF^ + d-yF^ + dpF^ = 0, (11) 
d a F af) = 0, (12) 
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In order to write in a concise manner the potentials of 
a continuous distribution of electrons in arbitrary mo- 
tion, one may assume that the customary equations, 
d a d a Ap = —Jp and d a A a = 0, hold. However, in 
this case Jp would represent the piecewise differentiable 
functions best approximating the 4-current of singulari- 
ties. The vanishing of the 4-divergence d a J a states the 
conservation of the net number of singular points. 

The repulsive Coulomb potential alone is not sufficient 
to get square integrable solutions. Tentatively, the sec- 
ond potential that is necessary for this purpose could 
be thought of as the quantum version of the vector po- 



tential of a point magnet: (A e /2) (<? x 



), where a 



represents the 4x4 Pauli spin matrices, and A e is 
a length-parameter measuring the strength of the cou- 
pling, not to be fixed before hand. The reader may wish 
to verify that the substitution of the matrix potential in 
equation is entirely equivalent to the substitution of 
the imaginary vector — i A e r~ 2 f. The covariant form of 
this 3-vector enters as part of the self-action 4-vector: 

S^(e) — -a Afj_ + i \ e cr^ ^-alo^ + i A e 9 M I 

-> -a&r-^-iXeir-^O) (13) 



The electromagnetic coupling constant in the equation 
above means that the electron was assumed to repel it- 
self with the same strength that two electrons repel each 
other. Now, the imaginary radial potential is also indis- 
pensable for the description of the electron-neutrino. As 
neutrinos may not have magnetic moment, perhaps this 
gauge invariant (a a p = d a ap — dp<J a — 0) new entity 
in physics is not related to magnetic moment but to the 
spin. For simplicity, the coupling constant A e was con- 
sidered to be the same for both particles. Eventually we 
will relate a to A e , which in turn could relate a to m v /m e 
if A e is written in terms of m e , m u , h and c. There are 
two options: 



(14) 



First solution: 

j + 1 — m 



n 1/2 



Vi 

'02 = - 

'03 = i 

4>A = i 



2(j + 1) 
j + 1 + m 



FY 



1.0 



j+2,m- 



(15) 



"I 1/2 



2(7 + 1) 

"I 1/2 



FY. 



\fy 1a 



j + m 



2j 
j - m 



GYj_^ m _i —>iGY 0fi 



-1 1/2 



2j 

Second solution: 



GY j _i m+ i — > 0. 



0i 

'04 



j + rn 
2j 

j - in 



1/2 



LY j _i^ m _i 



1/2 



2j 

j + 1 — m 



1/2 



20' + 1 

j + 1 + m 



KY J+hm _i 



i L Yo,q 


KY 



(16) 



1.0 



^/3 



-,1/2 



2(j + 1) 



KY j+hm+ x 



KY hl . 



The spherical harmonics are already normalized to one. 
The arrow indicates that the one and only value of the 
total angular momentum in conceptual connection with 
the problem is j — 1/2. (m = 1/2 is our option). The 
substitution of the sum of these solutions in equation 
QJ yields two identical systems of differential equations, 
each containing only one pair of radial functions, (F, G) 
or (K,L). We may leave functions (K, L) aside for a 
while and work with functions (F, G): 



d r F 



2F 

r 

d r G 



A e F 

KG 



E a 
he r 
E a 
he r 



G 



F 



(17a) 
(17b) 



where (3 = m u /m e . Interchanging the mass parameters 
would be the second option, but it does not work. 



III. THE ELECTRON SOLUTIONS 



The substitutions 



e 2 

E S > 



he m e 



F = F f exp(-a~ 1 /37s- 1 ), 
G = G f exp(-a- 1 /3 7 s- 1 ), 



E 



me* 



An electron in self-action is a state with definite 
energy which is taken into account with the factor 
exp(—iEt/h) in the solutions. Since the potentials are 
radial, the spatial part of the two independent solutions 
is to be written as shown Q . 



reduce the system to 

S - 2 d s (s 2 F^) = (l-s- 1 )aG^ (18a) 

dsG^ = -(l-r 1 )af t (18b) 

The reader may verify that the two independent 



4 



solutions, (F, G) and (f,g), can be written as: 



E 



E 



fc = 



a + > LrLI-sin- 

fc = l 



,.2fe 



E^(«' s ) 



oo 



E 



a 3 "" 1 (19) 

ft = o 

OO 

,-2 i ftf„\ „.2fc 



in.. 9j = 6 s +2^ fl(s)a A ", 

fe = fc = 1 



oc 



oo 



fc = 



,,2fc+l 



(20) 



Concerning the first solution, function G is the con- 
stant oo, and all the other functions as well as their first 
derivatives vanish at s = 1. Concerning the second so- 
lution, function /J = 6o s~ 2 , and all the other functions, 
as in the former case, vanish in a smooth manner at the 
classical electron radius. The procedure for the obtain- 
ment of solutions in powers of a 2 is iterative in nature 
and runs as follows. 

In the first solution, the first step is to insert G into 
(18a). Integration is the second step. The constant of 
integration serves to force the vanishing of /g at s = 1: 

s- 2 d s (s 2 F$) = (l-Oaao 
=> F$ =a(a Q /6){ S - 2 -3 + 2s). 

We proceed with the insertion of -Fjj into (18b). Inte- 
gration follows. The constant of integration forces the 
vanishing of G\ at s = 1: 

d s G\ = -(l-s- 1 )aF^ 
G\ = -a 2 (a /12)[s~ 2 - 2s" 1 + 6 In s + 9 - 10s + 2s 2 ] 

and so on. Notice that function G' cuts the s-axis only 
once at s 2 w a 2 /12. 

In the second solution, the iterative procedure starts 
with the substitution of /q = &o s ~ 2 into (18b), which 
yields function g^: 

d s g\ = -(l-s- 1 )ab s- 2 
=> gl = -a(b /2)(l-s) 2 s- 2 . 

The reader may follow the prescription to get functions 
fl and g\. With these functions at hand we may write 
down function g' G> s _1 which will be useful hereinafter: 

(g^G^s- 1 )s 2 = a^ + a 3 ^ + ... + a 2k+1 ^ k + ... (21) 

Co = a b [s _1 - 2 + s], 
fl = -Oo bo [s~ 3 - 4s~ 2 + 40 - 5s - 36s 2 
+4s 3 + 60s In s]/12 

Since functions: F, <?, gG and fF vanish in a smooth 
manner at s = 1, the physical properties of the electron 



should be represented with covariant densities made up 
of the product of the two independent solutions: take for 
example the 4-vector 



(22) 



where ip a and tpb are the two independent solutions. The 
time component of this 4-vector, 

(V>a7tV>6) - \\Y lfl \ 2 Ff + ^\Y 1 , 1 \ 2 Ff + Y 2 fl Gg (23) 

vanishes in a smooth manner at s = 1 ("ft is the unit 
matrix). However, we cannot consider the solution tp — 
for the region beyond s = 1 ; because doing so would make 
functions G and / look like strings with a loose end (see 
figure 1). 

Graphical discontinuities in the wave functions are not 
acceptable. This problem is solved by replacing the po- 
tential a Afj_ in favor of (m e cft _1 ) for the region be- 
yond the classical electron radius. The corresponding 
radial equations are 



d r F 



2F A P F 



d r G 



r- 
A f3 G 



E — m P c 



E 



G 



h c 



(24a) 
■ F (24b) 



Now we must set E = m e c 2 . With this consideration the 
external solutions, 

G = a exp(-A e r _1 ), F = 0, (25) 
f = b a r- 2 expi-Xer' 1 ), g = 0, (26) 

join with the internal solutions (19) and (20) at the clas- 
sical electron radius in a smooth manner. The external 
solutions have no further meaning since the correspond- 
ing bilinear products are null. 

The bilinear function (22) could be named interaction- 
density, the recipient of the action of the Coulomb po- 
tential. Accordingly, 



£ = 4 7r m e c 2 



\{s) (tpal^tpb) S 2 ds 



(27) 



would represent the electromagnetic inertia of the elec- 
tron. This definition, reminiscent of Maxwell's expres- 
sion (see equation (4)) is the formal substitution of the 
charge-density in favor of the interaction density. In 
Maxwell's theory the charge density determines the elec- 
trostatic potential, whereas in the fundamental theory 
here being proposed the interaction density is determined 
by two potentials endowed with independent existence. 



IV. FURTHER RESULTS 

The incorporation of the Dirac fifth matrix (75 in equa- 
tion (2)) into this approach enables one to write down the 
whole set of covariant densities H . 
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1 /So 


1 




(a) 


G 




/ \ -F 

/ \ 






(b) 




1 



\ (Ff+Gg) s" 1 
\ (c) 

v -^Ggs" 1 



FIG. 1: The graphs of the functions merely show their main 
features: all graphs vanish in a smooth manner at the origin 
of coordinates, (a) In the first independent solution, function 
F is (roughly speaking) a times function G. Function G cuts 
the s-axis only once at s 2 ~ a 2 /12. Function F vanishes in 
a smooth manner at s — 1. Beyond s = 1, function G is 
essentially constant, (b) In the second independent solution, 
function g is about a times function /; the former smoothly 
vanishes at s = 1. Beyond this radius, / decays essentially 
as 8~ 2 . (c) The products of the functions vanish in a smooth 
manner at s = 1. Beyond the classical electron radius the 
bilinear products are null. Thus, the energy of the electron is 
confined within its classical radius. 

The matrices, 

7t = diag [1, 1, 1, 1] , 

75 = diag [1,1,-1,-1], 

lx ly = diag 

Ixlylb = diag [i, -i, -i, i] . 



and the table below are particularly important in this 
regard. 



Pseudo Invariant 


1 


none 


Invariant 


1 




4-vector 


4 




pseudo 4-vector 


4 


(ipalxjy^b) 


anti-symmetric tensor 


6 


{^alxlylhlpb 



The first column shows the whole set of covariant 
densities. The second column indicates the number of 
components of each density. The third column shows 
the only component of the density that survives volume 
integration; the other components contain products of 
different spherical harmonics and their volume integral 
vanishes. 

The transformation properties of S z = {ipalxly^b) and 
M z = {^alxlyl^b) , suggest to interpret the former as 
the z-component of spin-density and the latter as the z- 
component of magnetic polarization. From the matrices 
above and equations (15) and (16) we can see that 5*2 is 
the same for the two independent solutions, (F; G, /, g) 
and (K, L,k,l). Though, M z , corresponding to the first 
solution is opposite to M z of the second solution. Since 
electrons and positrons with the same spin orientation 
have opposite magnetic moments and the energy of both 
solutions is positive, perhaps the negative energy solu- 
tions of the Dirac equation, as well as the philosophy 
built up around them, need to be re-examined. Curiously 
enough, this is not the only case where physicists have 
thought that space itself has strange properties: Special 
Relativity came to show the futility of the study of the 
hydrodynamic properties the ether that presumably per- 
vaded the whole universe. 

The graph of the density of the electromagnetic inertia 
(see panel (c) of figure 1) invites the volume integral of 
the function gGs^ 1 to vanish. This condition would in- 
terrelate a and (3 and the electron would look as a minute 
planetary system wherein two electro-quarks (named so 
because of the appearance of fractions 1/3 and 2/3 in 
(23)) retain all the electromagnetic inertia while revolv- 
ing around the singularity of the potentials. The covari- 
ant way of expressing this condition involves the Dirac 
fifth matrix 75: the density Ii = {ipal^b) is an in- 
variant and so is I2 = {ipalhi^bj ■ Therefore, the vanishing 
of the volume integral of Io(Ii — I2) is the condition we 
are talking about. 

For the purpose of evaluating this integral it is 
important to notice the following: 

1. If s ^> 2a _1 /3 = 7], function W = exp(— rys -1 ) is 
close to 1. 

2. Function W starts bending downward for values of 
s comparable to r\. 

3. For s <C i], function W is essentially null. 
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If r\ is sufficiently small we have: 

r 1 I 

/ W (1 + s + s~ 2 + s~ 3 )ds x 1 + - +r] 1 
Jo * 



Now, 



d n / Ws^ds 
Jo 

Therefore 

-l 



Ws^ds 



-T) 1 exp(— rf) w —77 x + l — 77/2+. 



In 77 + 77- ?7 2 /2! + ... + c (28) 







The constant of integration Co ~ —0.51 can be evaluated 
by plotting the graph of W s _1 (for 77 = 1) and com- 
paring the area under the graph with its analytical value 
given by the right hand side of (28) (evaluated for 77 = 1). 

When 77 is small enough, the dominant term of the in- 
tegral of W £1 (see equation (21)) is the first one. All 
these considerations imply that the vanishing of the vol- 
ume integral of G g s _1 is congruent with the following 
approximation 



-In 77 -0.51 -2 + 1/2 = (a 2 / 12 ) 7 ?" 



(29) 



Considering that m e — 511000 eV and that a — 1/137, 
we get that the mass of the e-neutrino is larger than 
1.75 eV but smaller than 1.78 eV. Hitherto experimental 
physics has set an upper limit of 2.2 eV for m„ |(|. 



V. THE ELECTRON-NEUTRINO SOLUTION 

Equation (1) provides an amazingly simple solution to 
the electron-neutrino if it is regarded as an electron de- 
prived of Coulomb potential. In this case the self-action 
would take this form: 



Sn(v) = -m v ch 1 u u + i X e <9 U I 



(30) 



We already know the corresponding radial equations, yet 
there is only one solution of interest: 



F 



s~ z exp(-^ s" 1 ), G = (31) 



The fact that it is not possible to cut off the solution at 
any particular radius suggests the Copenhagen interpre- 
tation of the square of the wave function: the probability 
of the electron-neutrino to interact with matter beyond 
one wave length from the singularity of its own potential 



u^ip^tp* s 2 dsj ij u fl ipj fl ip*s 2 ds 
= l-exp(-2/3 2 ) « 2(m„/m e ) 2 < lO" 10 



Neutrinos \i and r could be described with other 
coupling constants, (An and A T ). 



VI. THE PROTON 

The description of the proton can be carried out think- 
ing of it as a positron endowed with a Yukawa neutral 
meson field. Consider the scalar equation 

8 a d a ( = (m^cJT 1 )^, (32) 

its time independent solution 

( = r^ 1 exp[-(7n 7ro cft _1 )r] 

was used by Yukawa in his meson theory of the nuclear 
force [f| . In general, the nuclear force is attractive beyond 
the proton radius and repulsive for shorter distances, 
which prevents the nucleus from collapsing. Accordingly, 
the simplest self-action that we can think of is 



-n(iiu - a An + i \ e <Ju 



(33) 



where n is a dimensionless coupling constant which must 
vanish when the meson mass vanishes. The spin potential 
is assumed to be the same for the electron, the neutrino 
and the proton, all the offspring of neutron decay. 

The solution of equation (1) for this case can be 
obtained following procedures very similar to the one 
used for the electron case. Moreover, it is not nec- 
essary to use numerical analysis for the obtainment of 
a good approximation when n is one of the two sim- 
plest fundamental expressions: n = m Wo /m p ~ 1/7, or 



m. 



m p /(m 7ro + m p 



1/9. Actually, the correct 



value of m e /m p is obtained with the vanishing of the 
volume integral of G g (a Lj + n £) when n m 1/11. There- 
fore, it can be inferred that the proton self-energy was 
defined as 

rsa 

£ = 4irm p c 2 / (a + n ( u M ) {ip a 7u ipb) s 2 ds (34) 
Jo 

where sq is the radius where the right hand side of the 
corresponding radial equation vanishes. 

This model can be improved to get a better approxima- 
tion of m e /m p . However, we see a drawback of a different 
kind: the prevailing consensus is that protons are made 
up of three spin 1/2 subparticles. Perhaps an attempt 
to express the GUT's model in a covariant way should 
be made. How? assigning to each quark an individual 
bispinor wave function, yet all of them interrelated in a 
system of differential equations similar to (1). 

One more interesting aspect of scalar equations is that 
the substitutions = ip(i) exp[— i (mc 2 + E)th~ 1 ] 

and A/3 dp — c _1 4> dt in 



yield 



(d d / 3-2iaA d l3 )'il> = {m 2 c 2 h- 2 )^ (35) 



-{h 2 /2m) V 2 V + e 2 <t> [1 + (E/mc 2 )} ip 

-- Eijj[l + (E/2mc 2 )]. (36) 
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Neglecting E/mc 2 , which in atomic physics is much 
smaller than self-energy, the familiar Schrodinger equa- 
tion is obtained. However, the important point is that 
imaginary vector potentials are necessary to yield differ- 
ential equations with real coefficients. This is a feature 
of relativistic wave equations with the time-dependent 
factor exp(— i t ftr 1 ). In fact, we had this situation with 
equation (1) and the potential i\ e a^. 

As a final remark, the time-dependent factor 
cxp (— Ethr 1 ) in equation (1) is a logical step to start 
investigating spin 1/2 unstable particles. 
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